Dissipative Dynamics in a Quantum Register 
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A model for a quantum register dissipatively coupled with a bosonic thermal bath is studied. 
The register consists of A'^ qubits (i.e. spin | degrees of freedom), the bath is described by Ni, bosonic 
modes. The register-bath coupling is chosen in such a way that the total number of excitations is 
conserved. The Hilbert space splits allowing the study of the dynamics separately in each sector. 
Assuming that the couphng with the bath is the same for all qubits, the excitation sectors have 
a further decomposition according the irreducible representations of the su{2) spin algebra. The 
stability against environment-generated noise of the information encoded in a quantum state of the 
register depends on its su{2) symmetry content. At zero temperature we find that states belonging 
to the vacuum symmetry sector have for long time vanishing fidelity, whereas each lowest spin vector 
is decoupled from the bath and therefore is decoherence free. Numerical results are shown in the 
one-excitation space in the case qubit-dependent bath-system coupling. 
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I. INTRODUCTION 

The unavoidable interaction that each real-world sys- 
tem has with its environment is one of the major limita- 
tions to pratical realization of a quantum computer [Q. 
Indeed the outstanding potential capabilities of such a 
device rely heavily on the possibility of maintaining the 
quantum coherence in the system, and one of the typical 
effects of the coupling with the environment is to destroy 
phase relations between quantum states appearing in a 
linear superposition. This latter phenomenon is referred 
to as decoherence ^ : it can take place also when there is 
no system-environment energy exchange at all. To over- 
come this difficulty, in the last few years there has been a 
growing interest in the so-called error correction schemes 
P|, in which, by means of suitable encondings and mea- 
surement protocols, one is able to disentangle the system 
from the environment in order to recover uncorrupted in- 
formation. Another possible approach, pointed out in [Q, 
is to make use of symmetry to protect the information 
stored in a quantum state against environment-induced 
noise. In this case one, rather than to design states that 
can be easily corrected, looks for states that cannot easily 
be disturbed. The simplest system for which dynamical 
symmetry provides us these "safe" states is a collection 
of N qubits coupled all in the same way with the envi- 
ronment. Such a system can be thought of as quantum 
register S, analogous to those characteristic of classical 
computation. In this paper we study a model Hamilto- 
nian describing the exchange of elementary quanta be- 
tween the register and the environment, modelled by a 
bosonic bath B. The coupling with the bath is realized in 
term of the off-diagonal generators of a su{2) dynamical 
algebra The marginal dynamics of S is dissipative 
(i.e. the register energy is not conserved). The global 
Hilbert space decomposes in (dynamically) invariant sec- 



tors, characterized by their su{2) symmetry content as 
well as their number of excitations, as it will be spec- 
ified later. It should be emphasized that the focus of 
this paper is on the role played by dynamical-algebraic 
structures in providing collective states of the register in- 
trinsecally stable against environment-induced decoher- 
ence. The adopted physical model is, in a sense, generic 
and it is not aimed to describe a specific physical im- 
plementation of a Quantum Computer (as done instead, 
for example in |^) but a broad class of open quantum 
systems that could eventually turn out to be relevant for 
quantum data processing applications. 
In sect. II, after recalling the fundamentals of open quan- 
tum systems, the model is introduced and its general 
feature briefly discussed. In sect. Ill the associated 
Hilbert space structure is analyzed. In sect. IV the one- 
excitation subspace is studied, analytical as well as nu- 
merical results are presented for the fidelity and entropy. 
Sect. V contains some preliminary numerical results in 
the case of a qubit depending coupling with the bath. 
Section VI contains a number of conclusive remarks and 
perspectives 



II. THE MODEL 

Before introducing our model we begin by briefly re- 
calling a few basic facts about open quantum systems. 
Let Hs, Ti-b denote respectively the system and the en- 
vironment Hilbert spaces. We assume Tib to be much 
larger than Tig- The total Hilbert space is given by the 
tensor product Ti = T-is^ Tib- A state over Ha (a = s, b) 
is a hermitean non negative operator pa of End(7iQ.) with 
tr"(pQ) = 1. The manifold of the state over Ha will be 
denoted by Sa- The elements of Sa that are also projec- 
tors (p^ = p) provide the pure states. The set of 
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pure states generates Sa as its convex hull, furthermore 
it is in a one-to-one correspondence with TCa- Accord- 
ing to quantum mechanics, time evolution of the overall 
(closed) system is unitary, therefore if the initial state 
has the separable form p(0) — Ps Pb, {pa & Sa) then 
for any t > 0, the marginal (Liouvillian) evolution on Sg 
(open) is given by 

£t-Ss ^ Ss-.ps ^ tr*- (Ut p{0) U}), (1) 

where tr'' denotes the partial trace over Tib. The superop- 
erators {£t''}t>o are trace-preserving completely positive 
maps , that are the most general description of the evo- 
lution of an open quantum system. is not invariant 
under the action of {ff''}t>o, typically an initial pure 
state of the subsystem becomes mixed in a very short 
time scale depending on the strength of the interaction. 
This state can either eventually get pure again or not, but 
in any case an irreversible loss of the information stored 
in the initial preparation has occurred. It is important 
to notice that this mechanism is active even when there 
is no energy-exchange between the subsystems (i.e. the 
subystem Hamiltonians are constants of motion) at finite 
as well as at zero temperature. When an energy-exchange 
occurs we call the resulting dynamics dissipative. 
We introduce now the model. The Hamiltonian is given 
by iJ = Hs + Hb + Hi where 

N 

Hs=Y.'.al (2) 

i=l 

Hb = J2^kiblbk + l/2), 
fe=i 

N Nt 

Hi ^J2Y.(9krbia- +h.c.), 

i=i k=i 

here the erf's are the spin 1/2 Pauli operators (i is the 
qubit index), and the fo^'s bosonic operators. Hs (Hb) is 
the Hamiltonian of the register (bath), Hj the register- 
bath interaction. This model is closely related to the 
one known in the literature as the Dicke maser model 
1^. The latter, for generic iV, Nb, is not solvable and 
has a non-trivial ground-state phase diagram. In order 
to shed some light on the physics of this system we write 
down the equation of the motion for the Heisenberg op- 
erators 0{t) = U{ty O U{t). To simplify the expressions, 
it turns useful to perform the (unitary) transformation 
(7^ t-^ (7^ exp(±zeit), bk bkexp{—iujkt), whereby the 
Heisenberg equations then read 

^^^2j29k^it)bial (3) 

k 

^-Y.(9k^{t)bla- -h.c), 

j 



where gki(t) = gkieicp[i {tOk — £i)t]. By a formal inte- 
gration of the field equation and the substitution of the 
result into the spin equations one obtains 

' ^ ^ " ^ drK.At,r)a+ir)arit) - h.c. 

Here 

K,,{t,t')=J29k^{t)gl,{t') (5) 

k 

(bt{t)=2j29kr{t)bli0)at{t), 

k 

0.^w = E5H(OfcI(o)^rw- 

k 

This coupled system of non-linear integro-differential 
equations describes the dynamics of the spin subsystem 
in closed form. By means of the intermediation of the 
bath bosons each spin gets interacting with all the oth- 
ers it via a sort of time-retarded Heisenberg coupling. 
The information about the bath (dynamics as well as 
preparation) is contained in the kernels Kij, and in the 
operators {(f>f}. If = e, (j = 1, . . . , iV) and the bath- 
spin coupling is the same for all the spins, one has 

K,j{t,t') ^Y.lgkl'' exp[-t{e - ujk) {t - t')]. (6) 

k 

Under rather general assumptions this kernel is strongly 
peaked at t = t',U one has K,j{t, t') - 5{t - t') then (|) 
become a system of coupled non-linear differential equa- 
tions. Despite this strong simplification, also in this case 
the solution, due to non-linearity, remains difficult and 
one has to resort to numerical techniques. 
An alternative approach based on symmetry considera- 
tions, will be introduced in the next section. 



III. HILBERT SPACE STRUCTURE 

The Hilbert space is given by the tensor product Ti ~ 
'Hf^®nf^\ being Hs (Hb) the two (infinite) dimen- 
sional single spin (boson) space. The coupling of the 
spin system with the bosonic bath is described by the 
hamiltonian Hj is such that the raising (lowering) of one 
spin state is associated to the destruction (creation) of 
one boson. From this follows that the system admits the 
constant of motion 

N Nb 

1=1 k=i 
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The eigenvalues of X give the number of elementary 
(spin as well as bosonic) excitations over the reference 
state |0) = |0)s (g) |0)b. The latter is a lowest weight 
vector for the spin as well as for the boson algebra: 
|0) = ^fe |0) — 0, Va, fc. Its energy is set equal 
to zero. The Hilbert space splits into invariant eigen- 
spaces of I, 7i = ©/Ti/; an elementary combinatorial 
argument shows that the dimension of the /-excitations 
space Hi, (/ e N) is given by 



di = 



min{N .1) 

E 



I-l + Nb-1 
Nb-1 



(8) 



li N ^ Nt = 1, one has do = 1, d/ = 2, (/ > 1), the 
model reduces to the exactly solvable Jaynes-Cummings 
model of quantum optics ||9|. It is worth noticing that 
the general spin-boson model considered in the litera- 
ture on quantum dissipation, [usually addressed in the 
framework of the Feynmann- Vernon influence functional 
(see 1 10 for a review)], due to the presence of terms 



bk fJ^ , h], (Tj , (neglected here in view of the rotating wave 
approximation) does not conserve X, spoiling the associ- 
ated dynamical decomposition of the Hilbert space, on 
which our subsequent numerical analysis relies. Never- 
theless, since in this papers we are interested only in the 
role played by collective effects in stabilizing a quantum 
state, this restrictions does not result in any severe loss 
of generality. 



A basis for Hi is given by \ip)^ 



a, 



(g) |fc), where 



n<io)., \k)^iibi\o)t 



(9) 



Where n = 1, . . . ,niin{I, N), a E C{N,n), k E 
C'{Nb,N - n), C(ri, fc) {C'{n,k)) denoting the set of 
the combinations without (with) repetitions of n ob- 
jects k by k. Following the general scheme of Q we 
specialize hereafter the model assuming the parameters 
{9ki}, (i = l,---,-^), independent of the qubit 
replica index i . The first assumptions follows simply from 
the fact that qubits are replicas of the same system. The 
independence of the coupling constants on the qubit in- 
dex is analogous to the so-called Dicke limit of quantum 
optics it holds - for example - when the typical bath 
wave-lengths coupled with the register are much greater 
than the distances between the qubits. In this latter case 
the environment is no able to probe the internal structure 
of the register: as long as the dynamics is concerned it has 
an effective point-like topology. The common value e of 
the qubit "magnetic" fields {ei} will be chosen as the unit 
of the energy scale; analogously the "sound velocity" of 
the boson will be set equal to one, so that their dispersion 
relation reads LOk = k, {k = 2'Kn/Nb, n — 1, . . . , Nf,). 
The Hamiltonian can then be written as 



being S" = X^^i '^j'^ (o^ = -2,±) global spin operators, 
spanning a Lie Algebra sm(2), and B = '^f^gkbk- 
The fact that in (|l^) the spins appear only trough the 
5"'s means the all the qubits are treated symmetrically: 
the dynamics allows only for coherent excitations of the 
computational (spin) degrees of freedom. This is, from 
the algebraic point of view, a very strong constraint: the 
dynamics gets invariant under the action of the symme- 
tryc group Sn of the qubit index permutations. This pro- 
vides us one more constant of motion. Indeed from ( p^ ) 
immediately follows that the total spin 5^ is conserved 
and Hi splits according the SM(2)-irrep. The multiplic- 
ity of each irrep associated with the total spin quantum 
number S is given by 



n{S,N) 



Nl (2S'+1) 



{N/2 + S + iy.{N/2- sy: 



(11) 



One finds the following decomposition in invariant sub- 
spaces 

N/2 n{S,N) 

'^liS^r), (12) 

S=S™(JV,/) r=l 
min(I-N/2,S) 

niiS, r)= |/, S, r, S') (g> Hb{,Nb{I, S')), 

where Sm{N, I) — max {N/2 — I , s), { s — ior N even, 
and s = 1/2 otherwise), Hb{N) denotes the eigenspace, 
in Hb, of Nb ~ Hk corresponding to the eigenvalue 
N, iVh(/,S'^) = I- N/2 - S'^ |/, S, r, 5^) is a simulta- 
neous eigenvector of X, S'^ , associated respectively to 
the eigenvalues /, S {S + 1), . The reference state |0), 
belongs to the subspace (BiHi{N/2), with maximal total 
spin eigenvalue S = N/2 : this subspace will be denoted 
by TY^J"", and referred to as the symmetric subspace. If 
= H'y"' n Hi, one has 

™ = span{(5+)"|<-")) I n = 0, . . . ,/} (13) 

min{N,I) 

dim(7if' ") = J2 

1=0 



I-l + Nb-1 
Nb-1 



H ^eS'- +BS+ + S- 



Hb 



(10) 



The orthogonal complement of Hf"^ will be denoted 
by Hf. It is the direct sum of all the sectors with 
non-maximal S'^-eigenvalue. Before ending this section 
we notice that an additional term of the form H' = 
S^^^Wk{bk + b\), would destroy the u{l) symmetry 
generated by X, but not the su{2) structure. Such a 
term, considered in Q , does not correspond to an energy- 
exchange but is a source of pure decoherence. Since our 
analysis relies on the invariant decomposition H = ®iHi 
this term has been omitted here. 
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IV. 1-EXCITATION SPACE 

Due to the field-theoretic nature of our model, the di- 
mensionality formula clearly shows that for increas- 
ing excitation number / the problem of diagonalizing 
H becomes rapidly intractable. In particular, a finite- 
temperature analysis (arbitrary number of excitations) 
is very difficult. Nevertheless one of the interesting fea- 
tures of quantum noise is to be active also a T = 0, 
thanks to vacuum fluctuations. This latter issue can be 
addressed by exact numerical means, without an artifi- 
cious truncation of the bosonic space, by noticing that the 
one-excitation space T-6'^\ has dimension c?i — N + N^, 
that is only a linear function of the total number of de- 
grees of freedom. The basis {iV'i^jj)} is given by |a) = 

a+ |0),(a = 1,...,^) and |fc) = bl |0), (fc = 1, . . . , iV^). 
Equation ( p^ in this case reads 

N-l 

Hi =Hi(7V/2,l)0Hi(iV/2-l,r), (14) 

The symmetric space (S — N/2) is Ni, + 1-dimensional 
and it is spanned by the vector = N'^/'^S^ |0), 

and by the whole set The subspace corresponds 

to 5 = N/2 — 1. An orthonormal basis of Ti^ is given by 

N 

^ S+ |0), S+ EE N-'^'Y.^'''^^ (15) 

where, k = 2rm/N, n = 1, . . . ,N — 1. Since in the 
bosonic vacuum factorizes, this subspace, when neces- 
sary, will be identified with its projection over Tig. Now 
we observe that the vectors {|^fc)}, are annihilated by 
S~ , as they have minimal S^-projection, but also by the 
{bk}, as they have empty boson sector. From this follows 
that Hi \(j)k) = 0, (Vfc) therefore Tif is decoupled from the 
bath; it is an energy eigenspace with eigenvalue E = e. 
In terms of evolution superoperators, if p is a state over 
Til, we have ithe fixed-point relations f °(p) = p, (i > 0), 
where £^ denotes the superoperator associated with the 
bath- vacuum density matrix |0);,(0|{,. The states over TL^ 
are unaffected hy the decoherence induced hy couphng 
with the bath vacuum and can therefore to encode in- 
formation in a safe way. The space Ti^ is noiseless only 
at zero temperature; for finite temperature the |<^fe)'s get 
mixed with all the vectors belonging to the same su{2)- 
irrep, making the induced dynamics non unitary. It is 
important to notice that, for N 2, this states are not 
the noiseless ones introduced in [Q, as the latter are as- 
sociated with spin singlets (i.e. are annihilated by 
and S~^) and are decoherence-free at any temperature, 
whereas the |0fc)'s belong to N — 1-dimensional su{2) 
multiplets. The spectrum in the symmetric subspace can 
be obtained by resorting to exact numerical diagonaliza- 
tion of H , that provides the eigenvectors and eigenvalues 



On the other hand the spectrum in Hl^"^ is given by the 
A^b -I- 1 zeros of the expression [0 

I |2 

Pn,nSE) ^E-e-Nj2 (16) 
fc=i 

that corresponds to the analogous single spin problem 
with rescaled coupling ^/N g^- This follows from 

the symmetry constraint that makes \tp^y"^) the only 
state coupled with the bosonic modes. Let It/iq) = 
Si l*^*)' i^i — {4'i\'4'o)) be the initial state; at i > we 
can write, in terms of the chosen basis 

dl 

\m) ^ e-^*|^o) - E ^7(i) l7) e n, (17) 

7=1 

C,(i) = Ec°c;e-^-*, (c;^(7|0.)). 

1=1 

The marginal density matrix is given by Ps(t) — 
tT''\il^{t)) {ip{t)\. By using the relations 

trV)(a'l -^^|0)s(OUa-,, tr'\k){k'\ ^ Skk'\0)s{Q\s, (18) 
trVXfcl =tr''|yfc)(a| =0, 

one obtains 

N 

Ps{t)^ Ca{t)Cc.'{t)\a){a'\ (19) 

aa' — l 

+ \o)MsJ2\^kit)\'. 

fe=i 

The first (last) N (Nf,) terms describe a sector with a re- 
versed (excited) spin (boson). The marginal density ma- 
trix can be readily diagonalized, simply by observing that 
it can be written in the form ps (t) = Pi (t) \ips (t)) {ips (t) \ + 
Po(i)|0),(0|„ where 

1 ^ 

\Mt)) = -j== E ^"(*) I") e ^f"^' (20) 

V^lW a = l 

and Po{t) = 1 - Pi(i) = Ef=i |Cfc(i)l'- The von 
Neumann entropy of S is therefore given by Ss{t) = 
— Yld=o -fi(0 log2 Pi{t). A completely symmetric expres- 
sion (q fc. Pi <-> Po) is obtained for the bath marginal 
density matrix pt, — tr^(p), from which follows that 
Sb = Ss, furthermore we have S{S\B) = S{B\S) = -Ss 
for the conditional entropies, and S{B : S) = —2 Ss for 
the mutual entropy; this is a consequence of the purity 
of the overall system-bath state. 

In order to study the corruption of the information stored 
in a pure quantum state it appears useful to study the 
following quantity, called (input-output) hdehty 

Fit)^{ro\Psit)\ro)- (21) 
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The fidehty measures the overlap between the initial state 
\iPq) and the evolved one. In the following we shall be in- 
terested in the evaluation of (|2|) with initial data in Hi 
of the form \^o) = m \0)b, (p.(0) = |V'g)(V'ol)- For 
such initial preparation, if \tps{t)) — Pi{t) \ips{t)), one 
can write 

F{t) = {Mt)\Ps{0)\Mt)) (22) 

N 
a=l 

The " decoherence" function D{t) is related to the de- 
cay of the off-diagonal elements of Ps{t); indeed if, at 
i = 0, we prepare the system in the state 2~^/^ (|0)s + 
IV'o)) 10)6, where 5*^ IV-q) = (1 - N/2) IV'o), it is imme- 
diate to check that {ilj^\ps{t)\0)s = D{t). The mech- 
anism responsible for the energy exchange induces also a 
dephasing between the zero and one excited spin states: 
dissipation is associated with decoherence. 
In figure ([|) is shown the behaviour of F{t), for different 
values of the coupling with the bath, with initial data 
li^o) = IV'*^™)- lu this figure and in the subsequent ones 
is chosen as the time unit. These results are ob- 
tained by exact diagonalization of H in Tii, that provides 
the dynamical functions {C^(t)}. For strong bath-system 
coupfings F{t) develops osciUatory structures, due to the 
back and forth exchange of energy between the system 
and the bath. We report the simulations for weak cou- 
plings, since it is the case physically relevant. Further- 
more, since we are essentialy interested in the role played 
in the large times dynamics, by the symmetry structure 
of the initial state, and not in a detailed description of 
the bath-system coupling, we have choosen gk = go, (Vfc). 
From the point of view of the energy- information loss the 
latter choice is the worst case in that each qubit is coupled 
equally well with all the bath modes, which should not, 
of course, be the case in real systems. The fidelity in this 
range of coupling parameters and for intermediate times, 
vanishes in exponential way F{t) ~ exp(— t/r). The re- 
laxation time T, which turns out to be inversely propor- 
tional to -^X^fclfffcl^' time scale over which the 
dissipative process takes place. The real and immaginary 
parts of D(t) have an exponential damping modulated by 
oscillations over a time scale e~^. For very small times a 
naive perturbation up to the second order in Hi shows 

that indeed F(t) l~t'^/2 N A, where A = J2k 15^1^- 
Of course this process is nothing but the relaxation of the 
excited spin, whose energy is transferred to the environ- 
ment; for sufficiently large times one finds 

Ut\r''n'»\0)b^\0)b'E>\i'b), (23) 

where jV'b) = Cfc |fc), is a superposition of all one bo- 
son states. 

Some words of caution are now in order. The model un- 
der consideration is nothing but a multi-mode generaliza- 
tion of the Jaynes-Cummings model with many atoms. 



In analogy with the latter, for long time scale, t > tg, it 
exhibits a complex pattern of collapses and revivals . 
Furhermore since in each excitation space we have only 
a finite number of degrees of freeedom the phenomenon 
of the Poincare recurrences is also present for t > tn. 
In the following we will show results for t <ti tc, tn, in 
other terms we assume that, thanks to the great number 
of bosonic modes and the weak coupling, the physically 
relevant time-scales are much smaller than the ones at 
which this more complex behaviour appear. The energy 
exchange of the register with the bath can then be con- 
sidered irreversible. 
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FIG. 1. Fidelity as a function of time for jV'o) = |V''^'">- 
The coupling function is gk{i) = g Vk,i, {N = 2, Nt ^ 200). 
The time unit is . 



Suppose now that |'i/;g) = Cg IV'^^™) + Ca where 
|'0") & Ti-i, is a normafized vector, and |cap + \cs\^ = 1- 
Also, for the sake of concreteness, and without any loss 
of generality, let us consider the case N — 2. The vectors 
I'lp'^y"^), 1-0"), are given respectively by the two Bell-basis 

states \r^"^) = ^(1 n)+i it)), i^-^) = ^(1 n)-i it)). 

The initial marginal density matrix of S is therefore given 
by 

p.(o) = icp \r'^^){r^"^\ + |cj2 \r){r\ m 

By using our previous result for the symmetric initial 
state, and the fact that 1-0") is an energy eigenstate, it 
follows easily from (^3|), for t large enough, that 

£?{\r){r'"^\) = £?{\r''n{r\) = o, (25) 
£?i\r){r\) = \r){ri snir^mr'n) = io).(oi.. 

Therefore the large times density matrix is given by 

PsO,\Ca\^\r){r\ + \Csms{0\s, (26) 

from which straightforwardly follows for the fidelity the 
behaviour 
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F^\cX = (1 - \cs?? = (1 - |(V/^"|V'o)n'- (27) 



In other terms: the final state depends on the initial 
preparation symmetry content; a complete corruption of 
the initial information is obtained only if the initial state 
belongs to the vacuum S'^-sector Hi{N/2), so that the 
smaller is the projection over it, the closer to one is the 
fidelity. 

The extreme case is |V'o) € 'H^, F{t) = 1, in which 
there is no relaxation at all. In the intermediate situa- 
tions the spin system remains partially entangled with 
the environment and its state never gets pure. This sit- 
uation is illustrated in figures and (^) where fidelity 
and entropy are shown as functions of time in the case of 
\^^) = \a)^ for M = 1, 2, 3. In this case it 

is trivial to check, by using equations (p^) , and (|2^) , that 
F ~ (1 - M/Nf, and S, ~ {M/N - l)log2(l - M/N) - 
M/N\og2{M/N). Notice that if |csp = 1, one has the 
complete de-excitation of the spin system, therefore the 
initial state j^/;*^'") is maximally entangled and the final 
state |0)s, with zero mutual entanglement of the qubits. 
The system undergoes energy as well as information loss. 
On the other hand if |cap = 1, the final, and initial, state 
IV'o) is maximally entangled: energy and information are 
conserved. 
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FIG. 3. Entropy as a function of time for 

\ro) = M-^^^J2iLi (M = 1, 2, 3; iV = 4, N, = 200). 
gk{i) = 0.01 Vfc,i. The time unit is e ^. 
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FIG. 2. Fidelity as a function of time for 

\ro) = M-^^^J2tU 1«>. =1,2,3; N ^ 4, iV^ = 200). 
Qkii) = 0.01 yk,i. The time unit is e ^. 



In the temporal range in which the decay of fidelity 
is exponential, the following relation between relaxation 
times holds: t{cs) |csp = t(1). 



V. REPLICA DEPENDENT COUPLING 



In this section we present some numerical results in the 
case in which the system-bath coupling depends on the 
qubit replica. If the coupling functions {gk{i)} and/or 
the qubit energies {e;} depend on the qubit replica, the 
total spin operator S'^ is no longer a constant of the mo- 
tion. This situation is expected to be more realistic than 
the one previously assumed in that the latter amounts to 
have a bath with an infinite (i.e. very large) coherence 
length. In this case the decomposition ( |l3| ) is not invari- 
ant: the dynamics results in a non-trivial mixing of the 
su(2)-irreducible sectors Ti.i{S,r). In particular one has 
S~Hi ^ 0, therefore the vectors IV"") can decay. In other 
words the loss of the symmetry constraint allows the 
dissipation-decoherence induced by the bath to invade 
the whole Hilbert space. We choose gk{i) = go cos(fcz/^), 
where is a parameter related to the bath coherence 
length (so that for if = oo we recover the results of the 
previous sections) . In figure (^ are reported the plots of 
F{t) with initial condition in for different if's. 
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FIG. 4. Fidelity as a function of time for £ Hi , 

with C = 10 (1), C = 5 (2), C = 1 (3). The coupling function is 
gk{i) = go cos(fci/^). N = 2, Nb = 200, 30 = 0-01. The time 
unit is e~^. 



Figure (g) shows the behaviour of F{t) for small times 
with ^ = 1. Notably one observes that the initial con- 
dition |-0q) S Tif, exhibits a faster fidelity decay with 
respect to {i/jq) G TCl^™', for short times t < tc- For 
longer times, with obvious meaning of the notation, 
FA{t) > Fsym{t). The numerical simulations show in any 
case that Fa > Fsym, the bar denoting temporal average. 
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replicas of a two-level system (i.e. a 7V-qubits register). 
The register environment is described by a bosonic bath 
consisting of N}, modes, with Ni, 3> N. Its coupling with 
the register is realized by the exchange of elementary 
quanta of energy. The resulting dynamics of the register 
is dissipative; in the weak coupling regime, energy and 
information are irreversibly lost into the bath. Even 
though the model is non-trivial, exact analytical as well 
as numerical results can be obtained thanks to the ex- 
istence of a constant of the motion (excitation number) 
that leads to a decomposition of the total Hilbert space 
in dynamically independent sectors. Assuming that the 
environment couples in the same way with all the register 
qubits one has a further splitting of the sectors according 
to the irreducible representations of the spin su{2) alge- 
bra. Each su{2) lowest vector is decoupled from the bath 
vacuum fluctuations and therefore, at zero temperature, 
is decoherence free. The smallest subspace in which one 
can have non-trivial physics is the one-excitation sector 
Til. The dimension of Tii scales linearly in the total num- 
ber of degrees of freedom, therefore a thorough analysis 
of the dynamics in Tii, can be performed by means of ex- 
act numerical diagonalization of the model Hamiltonian. 
The temporal dependence of quantity of interest, such as 
fidelity and entropy have been studied. The asymptotic 
behaviour depends on the symmetry content of the initial 
state. Smaller is the projection of the initial state over 
the vacuum su(2)-sector Ji^v™-^ greater is the fidelity. 
In particular a complete energy-information loss occurs 
only when the state belongs to Ti*^"'. Some numerical 
results for bath-system coupling dependent on the qubit 
are also presented. In this more realistic situation the 
su(2)-structure is unstable: dissipation and decoherence 
affects the whole Hilbert space and then safe encondings 
no longer exists. Nevertheless our results shows that, on 
long time scales, the average fidelity of the previously 
noiseless states is still greater of the one of the other 
states. This suggest that the symmetry-based protection 
of quantum state suggested in Q| can be valuable in the 
general case. This last issue, along with the necessary 
finite temperature generalizations, worth further investi- 
gations. 



FIG. 5. Fidelity as function of time for [V-o) = IV'"'""), 
(dashed line) and |V'o) G T^i, (solid line). The coupling func- 
tion is gk{i) = go cos(fci). iV = 2, TV^ = 200, go = 0.01. The 
time unit is . 
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